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We present a formalized perturbation theory for Markovian open systems in the language of 
a generalized Schrieffer- Wolff (SW) transformation. A non-unitary rotation decouples the unper- 
turbed steady states from all fast degrees of freedom, in order to obtain an effective Liouvillian, 
that reproduces the exact low excitation spectrum of the system. The transformation is derived 
in a constructive way, yielding a perturbative expansion of the effective Liouville operator. The 
presented formalism realizes an adiabatic elimination of fast degrees of freedom to arbitrary orders 
in the perturbation. We exemplarily employ the SW formalism to two generic open systems and 
discuss general properties of the different orders of the perturbation. 



I. INTRODUCTION 

After more than a century of intensive research, many- 
body physics is an increasingly thriving field describing 
most of the phenomena appearing in nature. Its goal is 
to understand the macroscopic properties of large col- 
lections of interacting particles (typically of the order of 
10 23 ) from their microscopic laws of motion. In typi- 
cal situations the dynamics of the ensemble is governed 
by a Hamiltonian H whose complexity and vast dimen- 
sion impedes a direct solution. However, many complex 
quantum phenomena can be understood solely from the 
low-energy spectrum of H, such as quantum phase tran- 
sitions, topological insulation and superconductivity, just 
to name a few. Therefore a common strategy of many- 
body-physics is the derivation of a perturbative effective 
Hamiltonian H e g, which approximates the low-energy 
spectrum of H and reduces the complexity of the problem 
by integrating out the high-energy degrees of freedom. 
One of the most prominent examples of the success of 
this approach is the connection between the Kondo model 
and the low excitation spectrum of the Anderson model, 
which has been established in 1966 pp. It was achieved 
by a formalized version of (quasi-)degenerate perturba- 
tion theory [2H1] , which is nowadays known as Schrieffer- 
Wolff (SW) transformation and which paved the way for 
a deep understanding of these two distinguished models 
of condensed matter theory [5], The many analytical 
and numerical applications of this perturbative tool in 
contemporary physics are far too numerous to list here 
exhaustively (e.g., [fjlflO]). 

Due to the inevitable coupling of a quantum system 
to its environment, a paradigm shift could be observed 
in quantum physics in recent years, as the description of 
open systems moved into the focus of the field. Many 
seminal works in the context of, e.g., metrology in the 
presence of noise [TTJ [T2], dissipative quantum phase 
transitions }13[ 114) , as well as dissipation assisted quan- 
tum state preparation and quantum computation [T51 - 
[T7] . appeared over the past years. The situation for 
open systems in many respects parallels the considera- 
tions above. For Markovian environments, the system 
dynamics are described by a non-hermitian Liouville op- 



erator C. In many cases one is interested only in the 
low excitation spectrum of £, which describes the steady 
state behavior of the system and comprises in many situ- 
ations the relevant dynamics, since higher excitations are 
typically negligibly occupied during the system's evolu- 
tion. One prominent example constitutes the emerging 
field of dissipative phase transitions, which is intimately 
related to the low excitation spectrum of the Liouvillian 
[PEW20) . Also the widespread method of adiabatic elim- 
ination of fast evolving degrees of freedom - routinely 
employed in the field of quantum optics - corresponds 
to the derivation of a perturbative effective Liouvillian 
describing to low excitation dynamics of a system. More 
formalized perturbative tools that accomplish the goal of 
deriving effective dynamics for open systems to second or- 
der have been developed for specific scenarios [2TI 1221 155] . 
However, the available tools for open quantum systems 
are far less advanced than their Hamiltonian analogs. 

In this paper, we present a formalized perturbation 
theory in the language of a generalized SW transforma- 
tion, which adapts this formalism for Markovian open 
systems. We consider the most general case of a Li- 
ouvillian operator that features an internal hierarchy, 
i.e., it can be divided into a unperturbed part and a 
perturbation L = Cq + eV. A non-unitary similarity 
transformation on C dresses the zero eigenstates of Cq 
(i.e., the unperturbed steady states) with higher exci- 
tation eigenstates according to the perturbation eV and 
by construction decouples exactly the corresponding slow 
and fast space, respectively. The projection of the trans- 
formed Liouvillian onto this slow space (spanned by the 
dressed steady states of £ ) reproduces the exact low 
excitation spectrum and describes the system evolution 
in the vicinity of the steady state. In analogy to the 
unitary SW transformation this effective Liouvillian L e g 
can naturally be expanded in orders of the perturba- 
tion parameter e, yielding a systematic perturbative se- 
ries of the low excitation spectrum and in particular the 
steady state properties. We stress the point that in con- 
trast to previous perturbative approaches, our formalism 
works with minimal assumptions on the specific nature 
of Cq and produces perturbative results to arbitrary or- 
der. The procedure formalizes the usual perturbative 
approach and corresponds to an adiabatic elimination of 
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the fast evolving degrees of freedom. 

The phenomenon that the effective low-energy Hamil- 
tonian derived from integrating out high-energy degrees 
of freedom, often features a higher complexity than the 
original one, led in closed systems to the concept of per- 
turbation gadgets [23-25 . Along these lines, the idea of 
dissipative gadgets, i.e., the engineering of dissipation for 
quantum state preparation and protection has recently 
been proposed [TBJ [THl 127] . The presented SW formalism 
provides a natural tool for designing dissipation accord- 
ing to the desired steady state properties. 

The paper is structured as follows. In Section[ll]we de- 
rive the generalized SW transformation for open systems. 
We show, that in the new basis a subspace of slow dy- 
namics decouples exactly from all fast degrees of freedom 
and we derive an effective Liouvillian within this sub- 
space in a perturbative series [Eq. ( 27 )] . Subsequently, 



The in general non-orthogonal projectors 



in Section \LU\ we employ the formalism in two generic ex- 
amples, presenting two alternative strategies to evaluate 
the formal expressions for the effective Liouvillians and 
discussing general properties of the different orders of the 
perturbation. Finally, in Section |IV| we summarize the 
results and provide a brief outlook. 



Ct-.XaGV 

Q = l-P= £ \r t )(h\. 



(4) 
(5) 



define the subspaces corresponding to the spectral sets V 
and Q. These subspaces arc in the following referred to 
as the slow (defined by P) and fast (defined by Q) space, 
respectively, according to their evolution under the action 
of C . 

We use this partition of the left and right eigenbases 
to introduce a block structure for arbitrary superopera- 
tors A : L(H) — > L(H), where L(H) denotes the space 
of linear operators acting on the system's Hilbert space 
H (an example for such a superoperator constitutes the 
Liouville operator C itself), 



.4 



A p A~ 
A+ AQ 



(PAP PAQ\ 
\QAP QAQj 



(6) 



Further we introduce block diagonal and block off- 
diagonal operators 



II. FORMALISM 

We consider an open system whose evolution is gov- 
erned by a Markovian master equation. The correspond- 
ing Liouville operator can be partitioned in a zeroth order 
term Co and a perturbation V 



p = Cp = (C a + eV)p, 



(1) 



where e denotes the dimensionless perturbation parame- 
ter. £q is a linear operator on the vector space of C dxd 
matrices (d is the dimension of the system Hilbert space). 
We introduce the set of left and right eigenvectors for the 
non-hermitian operator Co- 

Co \n) = A 4 \ n ) , (2) 
(k\C = \i{k\, (3) 

which are chosen to be biorthonormal (h\rj) — 5jj and 
generically satisfy the completeness relation \ r i) {h\ = 
11. The eigenvalues A^ are in general complex. 

Since Co is the generator of a universal dynamical 
map (i.e., a contractive semigroup), its eigenvalues fulfill 
Re(Aj) < 0. The generated maps are trace preserving, 
which guarantees that the kernel of Co is at least one 
dimensional. We partition its spectrum in two subsets 
V = {\ a \\ a = 0} ^ {} and Q = {A,|A 4 ^ 0} (throughout 
the paper we will refer to eigenvalues from the two sets 
and the corresponding eigenvectors with greek and ara- 
bic indices, respectively). The spectral gap of the unper- 
turbed Liouville operator is denoted as A = min(IAJ). 
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By construction the unperturbed Liouville operator is 
block diagonal in this basis Co = Co = Cq, while the 
perturbation in general contains both block diagonal and 
off-diagonal terms V = V D + V° . In analogy to the 
Hamiltonian Schrieffer- Wolff Transformation our goal is 
to find a similarity transformation 



C^rL = U^CU, 
such that the two subspaces decouple 



L° = 0. 



(9) 



(10) 



Being similar [Eq. (|9j)], the transformed (L) and origi- 
nal (£) Liouvillian share the same spectrum. In the per- 
turbative limit A > 2e||V|| [1], the eigenvalues of the su- 
peroperator L c jf = PLP (referred to as effective Liouville 
operator) reproduce the exact low excitation spectrum of 
C. The master equation 



A — L c ffp, 



(11) 



thus describes accurately the steady state properties and 
low excitation dynamics, i.e., the system evolution in 
vicinity of the steady state. In addition, the Schricffer- 
Wolff transformation offers by construction a natural ex- 
pansion of the effective Liouvillian in the perturbation 
parameter e. 
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In the following we generalize the generic procedure to 
construct the transformation matrix U for hermitian ma- 
trices (for a review see, e.g., [SIS]) to the non-hermitian 
case. It can be shown that Eq. (JToJ) does not uniquely 
define the decoupling operator U. We will here only con- 
sider the so called 'canonical' choice U = e , where the 
generator S is imposed to be block off-diagonal S D = 0. 
Other choices of S D are possible, which then lead to dif- 
ferent perturbation theory formalisms, as outlined in [3] 
for the hermitian case. Depending on the specific prob- 
lem, alternative gauge choices for S D may prove to be 
advantageous. The discussion of the properties of the 
various formalisms represents an interesting subject for 
future studies. 

To simplify the formalism we introduce a compact no- 
tation where the commutation with an operator A is ex- 
pressed via the superoperator A defined via: 



AB = [B, A]. 



(12) 



Here, B denotes an arbitrary operator. This notation 
allows for an compact representation of the similarity 
transformation Eq. (|9| 

L = U- l CU = e- s Ce s 

= C + [C,S} + ^[[£,S},S} + ... (13) 

OO ^ 

^ i\ 

i=0 

We partition the latter superoperator into its odd and 
even powers 



cosh(S) + sinli(S). 



(14) 



The convenience of the canonical choice now becomes 
evident. While the odd operator sinh(£) changes block 
diagonal to off diagonal operators and vice versa, the 
even powers cosh(5) respects that structure. Therefore 
we can rewrite condition Eq. ( 10 1 



By construction the projection of the zero order Liouvil- 
lian into the fast space QCq = CqQ — Q£oQ has full 
rank and its inverse is well defined. For simplicity we de- 
note (QCoQ) 1 = Cq 1 . For block off-diagonal operators 
X = X° the resolvent operator fulfills 



Ko(X£ )=X, 



(18) 



as can be checked straightforwardly. Applying this su- 
peroperator to Eq. ( 16 1 gives the conditional equation 



for the generating matrix S 



S = -eK SV D - eK Scoth(S)V ( 



(19) 



Having derived a formal implicit expression for the 
transformation matrix S, which renders the Liouville op- 
erator block diagonal L° = 0, we now derive a compact 
expression for the diagonal blocks LP in terms of S, V 
and £q. As before, the block off-diagonal structure of S 
allows us to write L D as combination of even and odd 
powers of the superoperator S 



L = lJ> {e s cf 



(20) 



= cosh(5')£ 
= L D - 



sinh(S , )£ t 
d cosh(£) - l^o 



tanh(S") 



sinh(5)£ <: 



where in the second line we used Eq. (|15|). Using the basic 
trigonometric relation sinh(a;) 
tanh(x/2) we find 



(cos 



T(x) — l)/tanh(x) = 



L = C D + tanh(Sy2)£° 
= £ + e ( V D + tanh(5/2)V° 



(21) 



As discussed in all the above hyperbolic transforma- 
tions are well defined for infinitesimal transformation ma- 
trices S, which is guaranteed for appropriate perturba- 
tion parameters e. We denote the perturbative correction 
to £q as 



L° = (e s C)° = 

& sinh(5)£ D + cosh(S)C c 



(15) 



sinh(5) 6rD n ,,A, r0 



S 



-SC U +cosh(S , )£ t 







SC D = -Scoth(S*)£° 
& SC a = -SeV D - eScoth(S)V°, 



(16) 



where in the last step we used CP = Co + fV D and 
C° = eV°. Eq. (16) can be solved formally using re- 



solvent operator techniques. Since by construction the 
'slow space' defined by the set V contains only unper- 
turbed eigenvalues A Q = 0, it is £ P = P£ = and the 
resolvent operator takes the simple form 



K (A) 



QC^AP 



PAC^Q. 



(17) 



W = e [V u + tanh(Sy2)V c 



(22) 



Since by construction CqP = PCq = the effective Li- 
ouville operator in the slow space is given as 



L off = PLP = PWP, 



(23) 



and all dynamics in that space are at least of first order 
in the perturbation. 

Expanding S in orders of the perturbation parameter 
e in a Taylor series 



S 



71=0 



n q 



(24) 



and using Eq. ( 19 ) one can deduce a recursive equation 



for the S n . With these results we can directly construct 
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the perturbative correction W via Eq. ( 22 1 order by or- 
der. 

The first few Taylor matrices read 

So = 0, (25) 
S, = -K a V° = V-Cv 1 - C^V+, 

s-2 = n a v D s 1 = -ftoV^oV , 



The corresponding expansion for the perturbative correc- 
tion matrix W = J2^=o e ™W n can be found via Eq. (22 ) 



Wo 


= o, 




Wi 


= v D , 










m 







(26) 



of operator expectation values with equations of motions 
that close under the action of Co (in finite dimensional 
Hilbert spaces this is always the case) , the effective mas- 
ter equation can be readily evaluated, even if the ancilla 
system is high dimensional and its dynamics complicated. 
We show that the effective Liouville operator to second 
order is always of Lindblad form |26j , implying a Marko- 
vian evolution of the system. 

If the zeroth order Liuovillian Co is simple, it is ad- 
visable to explicitly calculate the matrices V and Cq . 
Given these matrices, arbitrary orders of the perturba- 
tion can readily be evaluated according to Eq. (27). In 
Section III B we consider an example recently studied in 
the context of superradiance in solid state systems [28] . 
which features simple zeroth order dynamics. We calcu- 
late explicitly the effective Liouvillian up to third order 
and show that the typically neglected third order has 
significant impact on the evolution of the system. 



A. General Ancilla Setting 



In [1] a formal expressions for the n — th order as well 
as a diagrammatic technique has been derived, which can 
be directly applied to the case of non-hcrmitian matrices. 
Straightforward evaluation of Eqs. ( [26] ) and subsequent 
projection onto V yields the first orders of the effective 
Liouville operator in the slow space [cf. Eq. (23)] 



■if 



2 --PVQC^QVP 
f = v- Cq X V q C^V^ 



Lf = PV D P = V p , 



(27) 



--{V p ,V-£ - 2 V+} + , 



where we employed the notation introduced in Eq. ([6| 
and {^4, B} + — AB + BA denotes the anti commutator. 
Note that C^ reproduces the well know second order 
result of adiabatic elimination in dissipative systems |21j . 



III. EXAMPLES 

In this section we will exemplarily employ the formal- 
ism developed above in two generic situations and present 
two alternative strategies to evaluate the expressions for 
the effective Liouvillians of Eqs. (27). 



First, in Section III A we consider the general setting 
of an ancilla system, which undergoes fast (in general 
dissipative) dynamics and is weakly coupled to a system. 
We adiabatically eliminate the ancilla to second order, 
employing the SW formalism. The coefficients of the ef- 
fective Liouvillian are expressed in terms of ancilla time 
correlation functions, which can readily be evaluated us- 
ing the quantum regression theorem. If one finds a set 



In the following we consider an example of how to ap- 
ply the formalism in a generic ancilla setting. A system is 
weakly coupled to an (unspecified) ancilla system, which 
undergoes fast (dissipative and/or coherent) dynamics. 
The Hilbert space of the total system is the product of 
the ancilla's and system's spaces V. = Ha <8 "Hs ■ We as- 
sume that the evolution of the ancilla is governed by fast 
dynamics given by Cq. Co contains an arbitrary combi- 
nation of Lindblad and Hamiltonian terms 



C-oX 



k 



j k (L kX Ll 



{L\L k , X } + ) 



i[H 0iX ], (28) 



where both the L^'s and Hq act only on the ancilla 
space. Let us for simplicity assume that Co features a 
unique steady state Co<J ss — such that the projector 
on the space of zeroth order steady states can be writ- 
ten in the simple form P\ = <J S s €5 Tr^(x) = a ss ® /i 
[29] . In the last step we introduced the reduced density 
matrix /i = Tr^(x). The weak coupling to the system 
is realized by the most general Hamiltonian interaction 
eVx = -itE^! A a ® S a ,x\- A a and S a are arbitrary 
hermitian ancilla and system operators, respectively. For 
notational convenience we will suppress the <S>-symbol in 
the following: A (g) S = AS. The full master equation 
thus reads 



£-X = £-oX + (Vx- 



(29) 



Note that the example we consider here corresponds to 
the often encountered situation of a bipartite system with 
separation of timescales. The ancilla evolution occurs on 
a timescale much faster than the system evolution. Thus 
we can consider the system's evolution under the condi- 
tion that the ancilla has settled to its steady state. The 
method presented above represents a formal approach to 
adiabatically eliminate the fast ancilla dynamics. 
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1. First Order 



We first evaluate expression © : 



The first order in the expansion Eq. (27) can readily 
be evaluated 

Lf X = PVP X = -iJ2 P ^ S ^ a ^ ( 3 °) 

a 

= P (A a a ss [S a , n] + [A a ,a ss ]S a fi) . 



The second term vanishes, since the trace over a commu- 
tator is zero: P[A a , a ss ]S a n = <r aa Tr A ([A a , a ss ])S a ^ = 
0. Thus we find the first order of the effective evolution 

Lf x = <JssLf[i = -ia^Y^KA^S^fi). (31) 



Since we are only interested in the evolution of the re- 
duced density matrix [i = TrA(x) we can trace out the 
ancilla degrees of freedom and find the first order correc- 
tion of the system's evolution 



Lf» = -iY^[{A a )S a , 



(32) 



Expectedly, to first order, the system experiences merely 
the effect of the semiclassical values of the ancilla opera- 
tors. 



2. Second Order 

The second order of the effective Liouville operator 
gives rise to more involved dynamics. We calculate the 
exact expressions and prove its Lindblad form for arbi- 
trary ancilla dynamics Cq. 

For the effective system evolution to second order we 
have to calculate the expression 

Tr A (Lf X ) = -Ty a {PVQC^QVPx). (33) 

In order to avoid the direct computation of Cq which 
may be impractical for large ancilla systems and for an- 
alytical purposes, we express the inverse via the Laplace 
transform Cq 1 = — d,Te c ° T , and we find 

/>oo 

Tr A (Lfx) = / At Ti A (PVQe CoT QVP\) (34) 
Jo 

/•OO 

= / drTi A [PV(l - P)e c ° T {l - P)VPx] 
Jo 

= dTTi A (PVe CoT VP X ) }® 
- jH ) drTv A {PVPVP X ), }© 
where we exploited the property Pe CaT = e CaT P = P. 



© 



dTTi A PVe L -° 



Car 



^ ' Aj Sj, cr ss fi 



(35) 



dr Y Tr A PV (e c " T A 3 a ss [Sj , fj] 

3 

+ e c ° T [A^a^nSj) 
= H) 2 Y { / drT ^ {Ai^Aj*™) [Si, [Sj,fj]} 

+ J dTTT A (A^lA^ass)) [Si^Sj] 

= -E { (/ [Si, [Sj,iA] 



In the last step, we defined the time correlation functions 
in the usual way (AiAj(r)) ss = Tr A (Aie^° T Aja as ) and 
([Ai^jfr)]).. = Tr A (A ie c ° T [A 3 ,a ss ]) . 

In the same fashion, @ can be readily evaluated to 
the formal expression 



2 = 



Y (l°° dT ) <^*>.-<4*>" [Si> iSj,(4] • (36) 



h3 



In general both formal expressions ® and @ are di- 
verging. However their sum 

Tr A (Lf X ) = ®+® (37) 
E{(/ driAAiAA^)^ [Si, [Sj, »]] 

J dT([AAi,AA jT }) ss } [S^nsA. 

represents a converging and meaningful expression. We 
defined AO t = 0{t) — (0) ss , for arbitrary ancilla op- 
erators O. © cancels the diverging parts in © and 
renders the integral over correlation functions finite, by 
subtracting the (infinite) steady state value. 

Next, we show that the second order derived above 



[Eq. ( 37 )] is always of Lindblad form, meaning that it 
generates a completely positive, trace preserving map. 
According to Eq. (37) the system evolution to second 



order is entirely determined by the matrix 



A= / dr{AAiAA jT 



IT .»„ J i dr(AAAA*} 

i,3 



(38) 



which can be written as a dyadic product of the vector 
AA = (A^4i,...,AA„) T . In fact, Eq. (37) can be rewrit- 



G 



ten in the more familiar form 

Tr A (Lfx) (39) 



— ^(■A-A^t^SiS-i,^ 



1,3 

Here it is evident, that the hermitian part of A is re- 
sponsible for the dissipative part of the evolution while 
the anti-hermitian part defines the coherent evolution. 
One readily checks that ^ J2i j(A — A^)ijSiSj defines a 
hermitian operator. On the other hand, in Appendix [A] 
we show the positivity of the coefficient matrix A + A^ > 
0, which guarantees Lindblad form of the dissipative term 
of Eq. (39). Thus, the evolution of the system after adi- 



abatic elimination of the ancilla is physical and up to 
second order Markovian. 

Aside from this general result, we now show that the 
coefficient matrix A can readily be calculated without 
evaluating the respective integrals explicitly, by using 
the quantum regression theorem [30j . Let us assume the 
equations of motion for the mean deviations of the ancilla 
operator set {^4 a } close under Co 



dt 



(AA t )=M(AA t 



(40) 



In finite dimensional system this can always be achieved 
by extending the set {^4 Q |a = 1, . . . , k} to a larger set 
{A Q |a = 1, . . . , n} (n > k) which forms an operator basis 
of the ancilla Hilbert space. 

Under these conditions the quantum regression theo- 
rem allows for a simple evaluation of the relevant time 
correlation functions 



— {AA t AA* 
at 

^ (AA t AA* 



M(AA t AA*\ 
■■ e Mt (AAAA* 



(41) 



All eigenvalues of the Bloch matrix M. have a strictly 
negative real part (and thus M is invertible), since Co 
generates a contractive semigroup with (by assumption) 
unique steady state. Therefore the latter equation can 
be readily integrated yielding a simple expression for the 
coefficient matrix 



/>oo 

A^ = dr(AA T AA*) ss = -M^iAAAA* 
Jo 



(42) 



The latter expression can be readily evaluated for a given 
system and uniquely defines the effective second order dy- 



namics of the system according to Master Eq. ( 39 1 . As 
shown in Appendix[A| independent of the nature and dy- 
namics of the ancilla system, the effective Master Eq. ( 39 ) 



is of Lindblad form and gives rise to a Markovian time 
evolution of the system. 

We emphasize that in many situations the size of the 
minimal set of operators that close under Co (defining the 



dimension of A4) will be much smaller than the dimen- 
sion of the Hilbert space. For illustration, consider the 
case where the ancilla system is constituted by a driven 
and damped spin-J and {A a |a = x,y,z} are the usual 
spin operators. In this case Eq. (40) represents optical 



Bloch equations and the matrix Ai is of dimension 3. 
Therefore, although the dimension of Co maybe large, 
the calculation of all coefficients of the effective second 
order dynamics reduce to a trivial low-dimensional ma- 
trix multiplication. 

We remark that a similar setting to the one presented 
here, has recently been examined in the context of dis- 
sipative quantum phase transition using the generalized 
SW technique [j"8] . 



B. Mediated Superradiance: Third Order 

In this Section we examine a specific example of an 
ancilla setting as discussed above. It recently has been 
studied in the context of superradiance from nuclear en- 
vironments of single photon emitters |28l . A radiatively 
decaying spin (e.g., a spin pumped electron spin in a 
quantum dot or nitrogen- vacancy (NV) center [3TJ I3"2"] ) 
is weakly hyperfine coupled to a large spin environment 
(e.g., nuclear spins of the host material). After each pho- 
ton emission the electron spin can escape from the dark 
state of the dissipation via exchange of an excitation with 
the nuclear spin environment. Supcrradiant features in 
the photon emission originate from a collective enhance- 
ment of the hyperfine flip-flop interaction for highly sym- 
metric nuclear states. 

In the following, we derive the effective evolution of the 
nuclear system after adiabatic elimination of the electron 
spin up to third order. In contrast to the previous exam- 
ple, where we expressed the effective Liouvillian in terms 
of integrated time correlation functions (which were eval- 
uated using the quantum regression theorem), we now 
calculate an explicit matrix representation of the pertur- 
bation operator V, in the biorthonormal eigenbasis of Co- 
Given this representation, all orders can readily be de- 
rived by simple matrix multiplication. We will find that 
the third order in the perturbation significantly improves 
the accuracy of the perturbative evolution. 

The model we consider is governed by the master equa- 
tion 



X= (Co + V) X , 



where 



Vp 



iu) \cr T a ,xj , 
'1 



if) 



(a+r +a~I+) +<7+a-I z ,x 



(43) 
(44) 

(45) 



where I a = f^cf are collective nuclear spin oper- 

ators, while erf and a a are individual nuclear and elec- 



7 







\n) 


\U) 


#1 


Ai = 


In) = 111) 


\h) = III) + Itt) 


#2 


A2 = —7/2 + iw 


k 2 > = lit) 


\h) = lit) 


#3 


A3 = —7/2 — iu> 


|ra> = in) 


lis) = Itl) 


#4 


A 4 = -7 


In) = Itt) - III) 


\h) = Itt) 



Table I: Eigenvalues and left and right eigenvectors of Co. We 
used the simplified notation = \i) ® \ j) =tj|)i which 
are the vector representation of the basis matrices \i) (j\ of 
the electronic space (cf. text). 



tronic spin-1/2 operators, respectively (a — ±,z). The 
individual hyperfine coupling constants gi are normalized 
to J^ff? = 1. 7 and u> denote the photon emission rate 
of the electron spin and the hyperfine detuning, respec- 
tively. The nuclear and electronic system are weakly hy- 
perfine coupled with gyN <C 7, This model describes 
the superradiant evolution of an NV center coupled to a 
nuclear spin environment, as discussed in [25] , 



After the assignment \ = J2i,j Xi,j N) 01 

zJi.i Xi,j I 



X 



\j) (where {\i} is an arbitrary basis of 



the matrix vector space), the superoperators of Eq. (43) 
can be written in matrix representation as 



£0 =7 



1 



a~ ® (ct + ) t - - (cr+cr- ® 1L + 11® {a + (j-) T ) 



(46) 



(<j + <j- ® 11- 11(g) (o- + O t ) , 



and 



V = -iff 



i (0-+/- + + a+a-p^j ® 11 (47) 



1L I I (<7+/- +0--/+) +0-+CT-/ 2 



where the superscript T denotes the matrix transpose. 
Since £q acts only on the electronic space, it can be 
straightforwardly diagonalized in the basis of left and 
right eigenvectors 



(48) 



The eigenvalues and the biorthonormal ((h\rj) = 5ij) 
left and right eigenvectors are given in Table ||J 

The representation of the perturbation in this basis 



V U = <Ji|V|rj), 
can readily be derived and is given as 



(49) 
(50) 





(v p 


v- 


V = 






v+ 












-10/2 (/"®11-1L® (/") T ) 


-»5/2 (/+®1L-1L®(/+) T ) 


-i.g (l z <g 


11- 


11® (7 Z ) T ) 


ig/2 U® (7+) T 


ig 11® (7 Z ) T 







81- 


f U® (7+) T ) 


— z.g/2 J" ® 11 





-ig I z ® 11 


z«?/2 (J- € 


9 1L + 


1L® (7-) T ) 





-ig/2 7" ® 11 


t S /2 U® (7+) T 


-i.g (l z (g 


11- 


1L® (,P) T ) 



Note that for simplicity we denote operators [e.g., Eq. ( 52 ), according to Eq. (27 1. In the following we calcu 



Eq. ( 49 )] and their representation in the Cq eigenbasis 
[e.g., Eq. (51 1] with the same symbols. The inverse of £0 
in the fast space is simply given as 



Co 1 = (QAQ) = 




(52) 





I/A3 
I/A4 

All orders of the perturbation can now readily be de- 



rived from products of the above matrices Eq. (511 and 



late explicitly the first three orders of the effective nuclear 
Liouville operator L cS . The first order Lf ff vanishes, since 
the perturbation vanishes in the slow space, V p = 0. The 
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second order yields 



T Off 

L 2 



(53) 



=7off 



{i + ) T -Ut®(i+r) T + i+r®i) 



- zw cff [i + r o n- 1® [i + r) T ] , 

where 7 ff = -2Re(.g 2 /A 2 ) = g 2 j/[(-//2) 2 +uj 2 } andw off = 
Im( 5 2 /A 2 ) = -gV[(7/2) 2 +^]. 

In the standard representation this corresponds to the 
second order master equation 



t eff 

/i = L 2 fi =7 cff 



1 



r^i + - -{i + r,n} + 



(54) 




pert. 2 nd & 3 rd order 



1 1.5 2 

time [7" 1 ] 



x 10 



which describes the collective decay of the nuclear spins 
at rate 7 c ff, responsible for a superradiant evolution as 
discussed in [28] . It is of Lindblad form and agrees with 
the result derived using standard adiabatic elimination 
techniques. 

After having derived the matrix representation of V 
[Eq. (51)] the SW formalism allows us to readily evalu- 



ate higher order corrections according to Eq. (27). This 



contrasts the situation for the standard techniques of adi- 
abatic elimination, where the derivation of higher order 
terms is extremely tedious. The third order of Eq. (27) 
yields the more involved expression (note that V p = 0) 



Q 



(55) 



-\T1 



*7T2 [ J ~ ® (I + I Z ) T -1®(I + I Z I~) 
4A 



4A2A4 



'4A§ 



[I+IT 



11 - rr 



4a 3 a 4 



[rr® {i+) T -r® (i+rf] , 



All t erm s involve contributions of the term oc a + cr~ I z in 
Eq. (45). It enters via in the intermediate process 
of Eq. (55]). In contrast, the second order was entirely 
independent of this term. 

In Fig. [T] we compare simulations of the exact evolu- 
tion according to Eq. ( [43] ) (solid) with the perturbative 
solution up to second (dotted) and third (dashed) order, 
respectively, for a system of N = 100 homogeneously 
coupled nuclear spins. The spins are initially fully polar- 
ized in z direction. We find that the third order effective 
Liouvillian has significant impact on the accuracy of the 
approximation. 

It can be shown that the Liouvillian Eq. ( 55 ) is her- 



miticity and trace preserving. Furthermore, the density 
matrix of the system remains positive throughout the 
evolution, although the third order Liouvillian it is not 
obviously of Lindblad form. In the following we show 



Figure 1: Comparison of the exact evolution according to 
Eq. ( 43 1 (dashed) with the approximate solution up to second 
(dotted) and third (solid) order [Eq. (pjj) & Eq. (55}]. The 



figure shows the radiation intensity emitted from the electron 
spin (oc d/dt(I z )) for a system of iV = 100 nuclear spins, ini- 
tially fully polarized in z direction. The emission intensity 
shows clearly the characteristic superradiant burst [28]. Pa- 
rameters are ui = 7/5 and y/N g = O.27. The second order 
effective evolution (L 1 ^) shows significant deviation from the 
exact dynamics. Addition of the third order L| , improves 
the approximation substantially. 



that in the limit to = one can recover Lindblad form 
by adding terms of higher order in the perturbation. In 
this case, all non-Lindblad effects are of an order that is 
deliberately neglected. 



For ui = 0, Eq. (55) reduces to the simple expression 



(in standard representation) 

Lf» = i^j eS [i-txi + ,i z } + i-^j eS [i + rr,fi]. (56) 

The first term of L| can be combined with the first 
term of L§ ff 



7eff / pi 



.9 ,r- 



7cff e 2 i 



2 7 
r I-fiI+ 



(57) 



e 2 t 



0[{ 9 hf 



where we used relation Eq. (13). The term O [(g/7) 2 ] is 
of fourth order in the perturbation and can consistently 
be neglected. The resulting master equation up to third 
order then has Lindblad form 



=7eff 



.■II 1 jv:l 

I* T~ 



(58) 



I~IJLl + e 



■ 9 



+ i^- lcS [I+I z I ,n\. 
47 

The excellent agreement of perturbative and exact 
solution displayed in Fig. [T] supports the expectation 
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that similar arguments hold in the general case lo 7^ 
[Eq. (55)] and effects due to the non-Lindblad form of 
the effective Liouvillian are of higher order in the pertur- 
bation. This assumption as well as the possibility to use 
the gauge invariance under the choice of S D , in order to 
ensure Lindblad form of every order of the perturbation, 
is subject to future work. 



IV. CONCLUSIONS 

We presented a generalized SW formalism, which 
adapts the successful perturbative tool of Hamiltonian 
quantum mechanics to the case of open quantum sys- 
tems, whose evolution is governed by a Liouville opera- 
tor. In analogy to the coherent case, we derive a trans- 
formation, that decouples subspaces of slow dynamics 
from fast evolving degrees of freedom in a perturbative 
series. In comparison with alternative schemes for adia- 
batic elimination [5TJ [33] , the advantages of the pre- 
sented method are twofold. First, minimal assumptions 
on the specific type of the zero order dynamics Cq have 
to be made. The subspace to be eliminated can be high 
dimensional and undergoing involved dynamics. Second, 
our approach in principle is an exact decoupling scheme, 
and can be applied to all orders in the perturbation. This 
property is of particular relevance for instance in numer- 
ical studies of low excitation spectra of the Liouville op- 
erator, e.g., in the context of dissipative phase transi- 
tions (T5HMI, and for error estimation in the context of 
quantum information processing [27j . The SW formalism 
provides a natural framework for the engineering of dis- 
sipative gadgets, e.g., in the context of state preparation 
and protection [T51 [TH1 [2Z] ■ 

We employed the SW formalism exemplarily to two 
model systems and presented different schemes to eval- 
uate the expressions for the effective Liouvillians. In a 
generic ancilla setting, we proved that the effective evo- 
lution of the weakly coupled system is up to second order 
Markovian, irrespective of the specific realization of the 
ancilla. In a second example we demonstrated that - in 
contrast to the standard schemes of adiabatic elimination 
- higher order corrections can readily be derived within 
the SW framework. In this model, the third order correc- 
tion plays a significant role in the perturbative dynamics. 

Further, we point out that the freedom of a gauge 
choice in the derivation of the transformation matrix 
could lead to a set of alternative perturbative approaches, 
which, in analogy to the standard SW transformation 
and depending on the specific problem, could prove to be 
advantageous under certain conditions. Potentially, this 
gauge freedom could also be used to ensure Lindblad form 
of the higher order effective Liouvillians. Lastly, we men- 
tion the numerous theoretical results H] (e.g., linked clus- 
ter theorem, additivity of effective Hamiltonian), which 
have been derived in the context of coherent SW transfor- 
mations as well as different variations of the SW method 
(e.g., continuous SW |34j ) . which may have open system 



analogs. These questions will be subject to future stud- 
ies. 
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Appendix A: Lindblad Form of Eq. (39) 



In the following, we prove that the dissipative part 
of the second order effective Liouvillian Eq. (39) is of 



Lindblad form. For this we have to show positivity of 
the respective coefficient matrix 

A + A f >0 (Al) 
^>v*(A + A^)v> 0, VueC" 

Expressing A as the integrated dyadic product 
[Eq. d38|] we write 



v*{A + A^)v (A2) 
dr ((v* AAAA*v) ss + (v*AA T AA*v) s 

dr ((aal) ss + (a T cr^) sa ) , 



where we introduced the ancilla operator a = v* AA. 

Since the expectation values are evaluated in the an- 
cilla's steady state, the two-time correlation functions are 
invariant under a total time translation t: 



(aal) ss = {cJ t <r\ +T ) t 



(A3) 
(A4) 



We exploit that property in symmetrizing Eq. ( A2 ) in 



the time arguments. First we 'average' Eq. (A2) over a 
total time translation 



v*(A + A*)v 



(A5) 



t roo 

dt I dr 



o Jo 



o 



(<°t<4fi-)«i + {a t+T a\) i 



1 



to poo 



dt 



' dt' ((<rt<ri),B + (ffo-j}") 



h Jo 

where the new variable t' = t + r has been introduced. 
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Basic integral transformations lead to the expression 



into two parts denned by the parameter x: J Q ° dt 



^ pto pto 

•-- dt dt'{a t a\,) 

l Jo Jo 



(A6) 



} 



/•to poo 

, dt / dt' {o t ol) ss 

*0 JO J to 
rto poo 

dt' dt{a t a\,) ss . 

to Jo J to 



} ® 
} © 



The first term of the latter equation ( © ) is positive 
since 



© = (RR f )ss > 0, 



(A7) 



with R = (l/V*o) / '° dt a t . 

We show that the remaining terms vanish in the limit 
to — > oo, proving the Lindblad form of Eq. (37). We 
estimate 



© 



■*o poo 

dt / 

• t a poo 

< — I dt dr\ (era. 

to Jo Jta-t 



1 

to Jo 
1 



dr{aal) ss 



(A8) 



where we reintroduced the time difference integration 
variable r = t' — t and used the time translation sym- 



metry Eq. (A3). Next the integration over dt is divided 



Jo 

as 



to-x 



dt + J ' °_ dt. The first term can be upper bounded 



to 



t -x 



dt 



< 



< 



dr \(aal) ss \ 



(A9) 



t -t 



Jo 

oo 



to—x 

dt 

dr I (cral) 



dr 



The second term can be estimated as 



1 

T 



t poo 

dt dr I [era. 

to—x J to — t 

2 pto P oo 

<— / dt 

to J t -x 
X 



T )ss I 



(A10) 



drUaal), 



to 



dr 



() 



Using the quantum regression theorem one shows that 
the time correlation function (aa\) ss decays exponen- 
tially. Choosing the parameter x = \fi both the right 
hand side of Eq. ( A9 1 and Eq. ( |A10 ) vanish in the limit 
<o — > oo. In an analogous estimation one shows the van- 
ishing of the remaining term © which proves the posi- 
tivity Eq. \kl\ and thus the Lindblad form of Eq. ([39]) . 
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